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The phase diagram of the multivalent manganites Lai_ x Ca x Mn03, in space of temperature and 
doping x, is a challenge for the theoretical physics. It is an important test for the model used to 
study these compounds and the method of calculation. To obtain theoretically this diagram for 
x < 0.5, we consider the two-band Double Exchange Model for manganites with added Jahn- Teller 
coupling and antiferromagnetic Heisenberg term. In order to calculate Curie and Neel temperatures 
we derive an effective Heisenberg model for a vector which describes the local orientation of the 
total magnetization of the system. The exchange constants of this model are different for different 
space directions and depend on the density of e g electrons, antiferromagnetic constants and the 
Jahn- Teller energy. To reproduce the well known phase transitions from A-type antiferromagnetism 
to ferromagnetism at low x and C-type antiferromagnetism to G-type antiferromagnetism at large 
x, we argue that the antiferromagnetic exchange constants should depend on the lattice direction. 
We show that ferromagnetic to A-type antiferromagnetic transition results from the Jahn- Teller 
distortion. Accounting adequately for the magnon-magnon interaction, Curie and Neel temperatures 
are calculated. The results are in very good agreement with the experiment and provide values for 
the model parameters, which best describe the behavior of the critical temperature for x < 0.5. 
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I. INTRODUCTION 

Manganites remain one of the most studied classes 
of materials in modern condensed matter physics, with 
many unanswered questions to keep us interested in 
them. At the same time, their diverse properties still 
hold the promise of many potential applications, and 
while low temperatures and high magnetic fields required 
may limit their availability in consumer electronics, they 
still provide excellent opportunities for advancing nan- 
otechnology. This is why in a series of papers 1 - - — during 
the past few years we have examined different aspects of 
the physics of manganites and constructed and tested a 
model which correctly describes their magnetic proper- 
ties, such as the observed phases and the Curie tempera- 
ture. In the present paper we will apply this model to a 
"real- world" compound, such as Lai_ x Ca x Mn03, and we 
will try to determine the values of the model parameters 
which best correspond to the experimental observations 
for the Curie temperature. 

Mixed-valence manganites with perovskite structure 
were first described by Jonker and Van Santeni in series 
of papers in the 1950s. These materials can be regarded 
as solid solutions between end members such as LaMn03 
and CaMnOa, leading to mixed- valence compounds such 
as Lai_ x Ca x Mn03. The general chemical formula for a 
perovskite compound is ABX3, where A and B are two 
cations of different sizes and X is an anion that bound 
to both. In the case of mixed-valence manganites, we 
have two different A type ions, usually alkaline metal 
or lanthanoid/rare-earth ions such as Ca 2+ , Sr 2+ , La 3+ . 
B type ions are the smaller Mn 3+ or Mn 4+ ions and X 
is oxygen. The ideal cubic-symmetry structure has the 
manganese cation in 6-fold coordination, surrounded by 



an octahedron of oxygen anions, and the A cation in 12- 
fold cuboctahedral coordination. The perovskite struc- 
ture is shown in Figure [1] for both end compounds. 

One of the most studied mixed-valence manganites is 
Lai_ x Ca x Mn03, which can be examined across the whole 
doping range < x < 1. While both of the end mem- 
bers LaMnC>3 and CaMn03 are antiferromagnetic and 
insulating, the resulting mixed-valence compound shows 
a variety of magnetic and transport properties depending 
on the value of the doping x. Wollan and Koehler— were 
the first to study the types of magnetic arrangements in 
the whole range of composition for Lai_ x Ca x Mn03 and 
to organize the results into a scheme of structures and 
structure transitions. At small hole densities, including 
x = 0, in addition to the already known pattern B, which 
corresponds to ferromagnetism, they observed pattern A. 
It corresponds to planes with ferromagnetic alignment 
of spins, but with antiferromagnetic coupling between 
those plains, and is known as A-type antiferromagnetism. 
Around x ~ 0.8, C-type pattern was observed, which cor- 
responds to ferromagnetic alignment of spins along the 
chains, but antiferromagnetic in the plains. When x ap- 
proaches one, almost all of the manganese ions are in 
the 4+ oxidation state and the dominant pattern is of 
G-type, with antiferromagnetism in all three directions. 
The corresponding types of magnetic arrangements are 
depicted in Figure [2] 

As a typical member of the transition metals family, 
Mn has a incomplete 3d shell filled with four and three 
electrons for Mn 3+ and Mn 4+ respectively. Due to the 
crystal field splitting effect, the degeneracy of the five 3d 
orbitals is lifted and they are grouped into one triplet t2 g 
and one doublet e g . The triplet has lower energy because 
of the space orientation of the corresponding orbitals in- 
side the octahedron of six oxygen ions, surrounding the 
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Figure 1. (Color online) Unit cell of a perovskite, shown for LaMnC>3 (left) and CaMnC>3 (right). 




Figure 2. (Color online) Types of magnetic ordering, as introduced by Wollan and Koehlep^. Each sphere represents Mn ion and 
the sign/color represents the orientation of the z-axis spin projection. A- type describes planes with ferromagnetic alignment 
of spins, but with antiferromagnetic coupling between them. B-type is ferromagnetism, C-type corresponds to ferromagnetic 
alignment of spins along the chains, but antiferromagnetic in the plains, and G-type is the nearest-neighbor antiferromagnetism. 



central manganese ion. The population of the t2 g elec- 
trons remains constant and the Hund rule enforces align- 
ment of the three t2 g spins into a state of maximum spin 
S = 3/2. Then, the t2 g sector can be replaced by a local- 
ized spin at each manganese ion, reducing the complexity 
of the original five orbital model. The electrons from the 
e g sector however can move from ion to ion, maintaining 
the projection of their spin, and are called mobile elec- 
trons. The only important interaction between the two 
sectors is the Hund coupling between localized t2 g spins 
and mobile e g electrons. 

Zener— was the first to construct a model which cor- 
rectly describes the properties of transition metals with 
incomplete d-shells, such as manganese. He introduced 
three principles, which he believed govern the interaction 
between the incomplete rf-shclls of neighboring atoms 
and later, as a demonstration of his ideas, he applied 
these principles to successfully explain the connection 
between conductivity and ferromagnetism in mangan- 
ites, observed by Jonker and Van Santeni. He inter- 
preted ferromagnetism as arising from the indirect cou- 
pling of incomplete d-shells via conducting electrons and 
also sketched a possible mechanism by which conduction 
electrons move between manganese ions. In this so-called 
"double-exchange" mechanism the transfer must occur 
through the oxygen ion between the two manganese ions, 
as a simultaneous "double exchange" of electrons. This 



process is a real charge transfer process and involves an 
overlap integral between the manganese 3d and oxygen's 
2p orbitals. Because of the strong Hund's coupling, the 
transfer-matrix element has finite value only when the 
core spins of the Mn ions arc aligned fcrromagnctically. 
Zener also showed that ferromagnetism would never oc- 
cur in the absence of conduction electrons or of some 
other indirect coupling. 

The ideas behind the Double Exchange Model were fur- 
ther developed by Anderson and HasegawaZ, who showed 
that electron transfer between neighboring Mn ions de- 
pends on the angle between their magnetic moments as 
t e g = t cos (0/2). The transfer probability varies from 
one for 6 = to zero for 6 = ir and the exchange energy 
is lower when the itinerant electrons spin is parallel to 
the total spin of the Mn cores. Starting from their re- 
sult, de Gennes^ introduced the so-called "spin-canted" 
state as a possible explanation of the coexistence of fer- 
romagnetic and antiferromagnetic features. Another pio- 
neer theoretical study was carried out by GoodenougbJi, 
regarding the charge, orbital and spin arrangements in 
the non-ferromagnetic regime of the phase diagram of 
Lai_ x Ca x MnC>3. It was based on the notions of "semi- 
covalent bond" and elastic energy considerations. 

Simple theories based on the Double Exchange Model 
however have several problems. They overestimate the 
Curie temperature of most manganitcs, cannot describe 
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the huge magnitude of the Colossal Magnctorcsistancc 
effect, underestimate the resistivity values in the param- 
agnetic phase by several orders and cannot account for 
the existence of charge/orbital ordering, phase separation 
scenario and strong lattice effects/anomalies seen experi- 
mentally. Searching for a more elaborate models, Millis^ 
was the first to argue that the physics of manganites is 
determined by the interplay between a strong electron- 
phonon interaction and the large Hund coupling. His 
ideas were later expanded by Roder— , Zangi£, Yunokii^ 
and others. 

The Jahn- Teller distortions lead to lifting the degen- 
eracy of the e g orbitals, with energy difference of about 
1 eV. The oxygens surrounding the manganese ion read- 
just their locations creating asymmetry between the dif- 
ferent directions. This effectively removes the degener- 
acy of the e g orbitals. The lifting of the degeneracy due 
to the orbital-lattice interaction is called Jahn- Teller ef- 
fect. The splitting between the two e g orbitals is big 
enough to justify the use of the simple one-band Double 
Exchange Model for manganites as an acceptable approx- 
imation. More elaborate models should however consider 
the electron-phonon interactions explicitly and include 
the corresponding terms in the Hamiltonian of the sys- 
tem. The importance of the crystal structure is demon- 
strated by the result that, without the Jahn- Teller dis- 
tortion (i.e. for a cubic cell), LaMn03 would be a ferro- 
magnetic metal rather than an A-type antiferromagnetic 
insulator— ~—. The Jahn- Teller effect responsible for the 
orthorhombic distortion also leads to the stabilization of 
the A-type antiferromagnetism. 

The double exchange model with Jahn- Teller coupling 
is a widely used model for manganites. The procedures 
followed to obtain the essential features of the model 
are different: numerical studie a 13 i 17 , Dynamical Mean- 
Field Theory (DMFT)^^, ab initio density-functional 
calculations^, and analytical calculation o 18 i 21 . In spite 
of the common conclusion that Jahn- Teller coupling sup- 
presses the ferromagnetic state, the results are quite dif- 
ferent and do not match the experimental results. For 
example the calculated Curie temperatures are two and 
even three times larger then the experimentally mea- 
sured. Because of that it is important to formulate theo- 
retical criteria for adequacy of the method of calculation. 
In our opinion the calculations should be in accordance 
with the Mermin- Wagner theorem^. It claims that at 
nonzero temperature, a ID or 2D isotropic spin-S Heisen- 
berg model with finite-range exchange interaction can be 
neither ferromagnetic nor antiferromagnetic. We employ 
a technique of calculation^, which captures the essentials 
of the magnon fluctuations in the theory, and for 2D 
systems one obtains zero Curie temperature, in accor- 
dance with Mermin- Wagner theorem. The physics of the 
ferromagnetic manganites near the Curie temperature is 
dominated by the magnon fluctuations and it is impor- 
tant to account for them in the best way. In contrast 
with other theories, our model includes the contribution 
of both the localized and the mobile electrons, thus bet- 



ter accounting for the magnon fluctuations. With our 
previous papers^ successfully explaining the qualitative 
picture, we now turn to comparing quantitative results 
and determining the model parameters which best de- 
scribe the experimental data. 

II. EFFECTIVE MODEL 

To construct the effective model, we start with the two- 
band Double Exchange Model with added Jahn- Teller 
distortion and antiferromagnetic (Hcisenbcrg) term. The 
Hamiltonian reads 

H = H^e + Hjt + Haf- (1) 

The first term describes the hopping of e g electrons and 
the Hund interaction between the spin Sj of the e g elec- 
tron and the localized t2 g spin 

Hue = — ^ if/'cJ Q Cj +a /' Q — 2 Jh ^ s.; • Sj (2) 

iall'a i 

where Cy and cu a are creation and annihilation opera- 
tors for e g electron with spin a on orbitals d x 2_ y 2 (I = a) 
and d 3z _ r 2(l — b) at site i. The sums are over all sites 
of a three-dimensional cubic lattice, and a is the vector 
connecting nearest-neighbor sites. For the cubic lattice, 
the hopping amplitudes between the orbitals along the 
x, y, z directions are: 

^aa = — V^3t„(, = —v3t ba = 3t bb = t 

tla = V3t% = V3t* a = 3t& = t (3) 
t Z aa = t Z ab = t\ a = 0, il b = 4t/3 

The second term in ([2]) is the Hund interaction between 
the spin Sj of the e g electron and the localized ti g spin 
Si with 

s i = \ ^2 c ti a <f) c iiPi ( 4 ) 

la/3 

where a x ,a y ,a z are the Pauli matrices, and the Hund's 
constant Jh is positive. 

The Hjt part models the coupling of the e g electrons 
to the phonons: 

#jt - gJ2 (Q^i T xi + Q*r«) + ^ (Ql + Ql) (5) 

i i 

where 

T xi = ( C taa C iba + C^ Q Ci aQ ) 

a 

T zi = ( c iaa C ia.a ~ c tba C iba) (6) 

a 

are the so-called pseudo-spin operators, g is the electron- 
phonon coupling constant, and Q2i and Q^i are the Jahn- 
Tcllcr phonon modes. The second term in -ffjx is the 
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general quadratic potential for distortions with constant 
k. The important energy scale of the phonon-electron in- 
teraction is the static Jahn- Teller energy Ejt = 5 2 /(2fc)- 
To explain the experimentally observed antiferro- 
magnctic phases in manganites we will also need a 
Heisenberg-like antiferromagnetic term Haf, which can 
be represented as: 



(7) 



the Bose fields at and af are the true magnons of the 
system. 

An important advantage of working with singlet 
fermions is the fact that in terms of these spin-singlet 
fields the spin-fermion interaction is in a diagonal form, 
the spin variables (magnons) are removed, and one ac- 
counts for it exactly: 



(16) 



where we have different exchange constants Jaf for dif- 
ferent lattice directions. 

We switch to Schwinger-boson representation for the 
localized spin operators 

Si = \'52vt a (7 aP<PiP> with ^2ipf a <p ia = 2s. (8) 

a/3 a 

By means of the Schwingcr-bosons we introduce spin- 
singlet Fermi fields 



1 



'2s 



[ip+(T) CU1 (T) + tp+(T)ca2(T)] (9) 
[<Ph(t)cu2(t) - tpi2(T)Cill(T)] (10) 



and write the spin of the e g electron and the total spin 
of the system 



S* ot = S, 4 
in terms of the singlet fermions 



gtot 



(ii) 



Si, (12) 



where the above two formulas account for the fact that 
the spins of the e g and tig electrons are parallel. 

If we average the total spin of the system in the sub- 
space of the singlet fermions A and B, the vector 



(s! ot )/ 



(13) 



Invoking (|9T TTU)) we can also rewrite the JT term as 

HJT = 5 £ [Q2i + Kb***) (I 7 ) 

i 

+ Q* (*ta*^« - Kb***)] + \ E (Qli + Qld 

i 

The theory is quadratic with respect to the spin-singlet 
fermions and one can integrate them out to obtain the 
free energy of fermions as a function of the magnons' 
fields. We expand the free energy in the ferromagnetic 
regime in powers of magnons' fields and keep only the 
first two terms. The first term J-/0, which does not de- 
pend on the magnons fields, is the free energy of Fermions 
with spins of localized t2 g electrons treated classically. 
We fix the model parameters and consider this term as a 
function of the Jahn- Teller distortion modes independent 
on the lattice site s. One ca n then show that this term de- 
pends only on sjQ\ + Q\. If we represent Q2 — Q cos 7 
and Q3 = Qsin7, this allows us to fix 7 = (Q3 = 0). 
The physical value of the Jahn- Teller distortion is the 
value at which Tf§ has a minimum. In this way we ob- 
tain the distortion as a function of the density of e g elec- 
trons for different values of the Jahn- Teller energy and 
fixed Hund's coupling^. 

The Hamiltonian corresponding to the free fermion 
part has the form 

H { = V \e\ + e B , 

t / , I c qfc * ak ak'^ak^ak ak 
k 



identifies the local orientation of the total magnetization. 
Because of the fact that t2g-electron spin is parallel with 
e 9 -electron spin we obtain Mj 



fSi with 



M = S+\Y,({*t + *ti-*u + *fi)) 



(14) 



If we use Holstein-Primakoff representation for the vec- 
tors Mj with M as an "effective spin" of the system 
(M 2 = M 2 ), 



mKk^k - <mtc<k - <mtuKk , m 



where 



— 2t (cos k x + cos k y ) — s Jh — /U 
— 2t (cos k x + cos k y ) + s Jh — \i 



-bk 



1 



1 



2t — cos k x + — cos k y + — cos k z — s Jh — A 4 



Mj = Mji + iM j2 = ^J2M - ajaj a, 
M 



M n - iM, 



Mf =2M-afaj, 



afJ2M-af aj (15) 



( i cos k x + i cos k y + ^ cos k z J 4 

e(fc) = 2 ( _ 73 cos ^ + 73 cosfc ^^) (19) 
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are the dispersions of the spin-singlet fermions. In order 
to diagonalize this Hamiltonian we use a Bogolyubov-like 
transformation 



form in terms of the quasiparticlcs 

Hf = E [ E akfak A fak + E akfak 'fa 





<fl + 


vifbi 


*& = 


-v A f A 
u k Jak 


+ <fbi 
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ak 


<m + 


- B kfb B k 


^bk — 


_„B fB 
u k Jak 


+ <f B 



rpA f +A f A , rpB r+B f B 
"^bkJbk Jbk+^bkJbk Jbk 



(22) 



(20) 



where the dispersions of the newly introduced fermions 
/, /+ are given by 



with coefficients 



u k 



r R 



{l+x%)/2 R = A 1 B 

sign(e(fc)) v /(l-a;f)/2 



(21) 
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b ak 
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+ h 
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- e B ^ 2 



(23) 



The Hamiltonian can then be rewritten in diagonal 
I 



The second term in the Fermion free energy gives the 
spin-fermion interaction 



= - E 



<ij> 



L aa 

2s 



^+ B ^ B ) 

a j at J 



-2a) {^JK + ^ifO 



E 



<ij> 



<ij> 
bb 

2s 



^ iCT - 2s) (*+ A *£ + *+ B * 

E ^ [ (vtrVia - 2-) + ffj^ + *5 B «S + 



(24) 



<ij> 



2s) ( + + *+ A *A + 

/ l ai ^ bj ' ^ bi a] ' bj ai ~ * aj ^ OJ J 
bi y V aj x to r ^ bj i * at ^ ai x bj bj aj J 



The spin-fermion Hamiltonian is quadratic with re- 
spect to the magnons' fields and defines the effective 
magnon Hamiltonian in Gaussian approximation: 



H ^ = - E jaM « • M *+* 



(26) 



^eff = E P & ( a i ai + a i+a«i+a - <4 «i+a - Ana^) 

(25) 

Based on the rotational symmetry, one can supple- 
ment this Hamiltonian up to an effective Hciscnbcrg-likc 
Hamiltonian, written in terms of the vectors Mj 



where J a are the effective couplings 

" M 2 Af ' 



(27) 



The effective exchange constants depend on the lattice 
direction and are a sum of two terms. The first gives the 
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contribution of the antifcrromagnctic Hamiltonian ([7|l. 
rewritten in terms of the vectors M,-. The second term 
gives the contribution from the spin-fcrmion interaction, 
where p a are the spin-stiffness constants (j2"5]l . which also 
depend on the space directions a. They are calculated at 



where n^, are the occupation numbers for the quasi- 
particlcs respectively. 

The ferromagnetic phase is stable if all effective ex- 
change coupling constants are positive J a > 0. If one 
of them is negative, for example J v < 0, and the other 
two arc positive J x > 0, J z > 0, the stable state is 
A-type antiferromagnetic phase which has planes (a;, z) 
that are ferromagnetic (parallel moments), with antifer- 
romagnetic (antiparallcl) arrangement of the magnetic 
moments between them (see fig-O- 

The spin-stiffness constants, function of the dop- 
ing x, are depicted in Figure [3] for the pure Double Ex- 
change Model (Jaf = 0> -^JT = 0), across a broad range 
of values for the Hund's coupling Ju/t. In such sim- 
ple model we account only for the contribution of itin- 
erant electrons, and as a consequence p x and p y remain 
equal for all carrier densities. One can see that by in- 
creasing the value of Jn/t the spin-stiffness constants 
also increase, with p z being the dominant one. Con- 
sequently, the range of the ferromagnetic region, deter- 
mined in this case by the condition p a > for all a, also 
increases with increase of Jn/t. From the form of the 
spin-stiffness curves we can estimate a limit for the value 
of the Hund's coupling, needed to reproduce the observed 
phases. Since the ferromagnetic to A-type antiferromag- 
netic transition occurs around x = 0.08, it is clear that 
values of Jn/t < 10 are not acceptable. The minimum 
value for Jn/t that extends the ferromagnetic phase to 
and beyond x = 0.08 is around J-n/t = 12, in accordance 



zero temperature for fixed Hund's coupling, Jahn- Teller 
energy, charge density, and Jahn- Teller distortion deter- 
mined above using the technique introduced in previous 
papers^, and have the form 



(28) 

I 

with previous results 1 - ' 23 ' 24 . For a realistic models, val- 
ues as big as Jn/t = 20 should be used to ensure that 
we will be well above the limit, since the addition of 
both the antiferromagnetic and Jahn- Teller terms leads 
to suppression of the magnon fluctuations^ 3 -. Another 
important observation we can make about Figure [3] is 
that the maximum of the spin-stiffness curves is shifted 
to smaller values of x with increase of Jn/t. This in turn 
determines the behavior of the Curie temperature curves 
and will play important role in determining the model 
parameters later. We will come back to this observation 
once we have introduced the method of calculating the 
Curie temperatures. 

The simple model used to construct Figure [3] however 
cannot explain all the observed phases. For example, it 
is a well known fact that the CaMn03 end member of 
the Lai_ x Ca x Mn03 family is a G-type antiferromagnet 
(see Fig. [2]). In addition to the G-type antiferromagnetic 
phase, with decreasing the doping C-type antiferromag- 
netic arrangement is observed before reaching the more 
complicated CE phased— . Both the G-type and the 
C-typc antifcrromagnctic arrangements require the Dou- 
ble Exchange Model Hamiltonian ^ to be supplemented 
with antiferromagnetic term ([?])• However, from Figure 3 
it is evident that antifcrromagnetism cannot be added in 
a trivial way. If we have equal value of Jaf for all three 
lattice directions, to have G-type antiferromagnctism this 
value has to be larger than p z /M . As a result the other 
two spin-stiffness constants will be negative for all doping 
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values. Thus, if we want to describe the correct sequence 
of phase transitions A — >• B — >• CE — >• C — >• G phase, we 
need to have different values of Jaf in different lattice 
directions, with Jj[ F being the largest. To explain the 
existence of the A-type antifcrromanctic phase, having 
in mind that p x = p y for all doping values, one might 
also consider having different values for J^ F and J^ F . 
However, the addition of J arm- Teller distortions splits p x 
and p y , as discussed below, and we can describe the ex- 
istence of the A-type antiferromagnctic phase even when 
Jaf = ^af- For simplicity we choose to work with equal 
values for them, which in turn leads to equal values of the 
effective couplings J x and J v across wide range of doping 
values, up to the appearance of J aim- Teller effects. The 
importance of different values of Jaf in different lattice 
directions will become more evident when we discuss the 
phase portraits. 

To help illustrate the impact of the Jahn- Teller distor- 
tions on the spin stiffness constants, we have depicted 



them in Figure @] for the case of non-zero Ejt- In accor- 
dance with the experimental^ and theoretical^ results, 
we have chosen a value for the Jahn- Teller energy with 
diminishing contribution beyond x = O.f. One can see 
that the appearance of distortions is accompanied with 
a change of the slopes of the spin-stiffness curves. The 
distortion splits the p v (dashed) and p x (dotted) lines, 
causing p v to decrease rapidly, while p x is stabilized. At 
the critical density x = 0.08, p y becomes equal to zero 
and the system undergoes a transition from ferromag- 
netic phase to A-type antiferromagnctic phase. Thus our 
model requires the existence of Jahn- Teller distortions 
to correctly describe the transition from ferromagnetic 
to A-type antiferromagnetic phase. The other remain- 
ing spin-stiffness constant, p Zl also starts to decrease. 
Since the spin stiffness constants are a measure for the 
magnon fluctuations in the ferromagnetic phase, which 
in turn determine the Curie temperature, the addition 
of Jahn- Teller distortions leads to decreasing the Curie 
temperature in the ferromagnetic regime^. 




III. CRITICAL TEMPERATURE IN THE 
FERROMAGNETIC REGIME 



To calculate the Curie temperature Tq we utilize 
Schwingcr-bosons mean-field theory^. The advantage 
of this method of calculation is that for 2D systems one 
obtains zero Curie temperature, in accordance with the 
Mermin- Wagner theorem^. We start by representing the 
vector Mi by means of Schwinger bosons (cj>i a , 4>f a ) 



Figure 4. (Color online) Spin-stiffness constants p a as a func- 
tion of hole doping x in the presence of Jahn- Teller distortion. 



Mr = ^6+a 



2 2^ fi^o^W 

a/3 



= 2M 



(29) 
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Next we use the identity 



(30) 

and rewrite the effective Hamiltonian in the form 



where the constant term is dropped. To ensure the 
Schwinger boson constraint we introduce a parameter (A) 
and add a new term to the effective Hamiltonian pip . 

ffeff = H cB + A {<t>t &r - 2M) (32) 

We treat the four-boson interaction within Hartrec- 
Fock approximation. The Hartree-Fock hamiltonian 
which corresponds to the effective hamiltonian reads 



-^HF = g-E [ Jx (vtrW+e**) {ft+e^a' -Via' ) + Jy (Piafi+eyv) \<Pt+ey + J * (<Pi<T<Pi+e*<r) ( < Pi+e,a><P<<r') 

i 

-Jx {Via'Vi+e**') Vi+evvPio ~ J * {vt+e^'V™') vtrfi+W ~ J y (ftr'Vi+ey*') vt+e y a^<? 
- J V \<Pi+e v < r ' ( Pi<r') Pt^i+eyCT ~ 'h (vt'Vi+w') Vi+e.oVio ~ J z (<Pi+e, fta^+e z a + A^ {ft a fta ~ 2M) (33) 



It can be rewritten in more compact form as 

ia i 



that purpose it is convenient to represent it in the form 
A = ( u aJ a + M u a)- In terms of the new parameter the 



-boson dispersion is 



£ k = ^ [ u * ja (1 - COS fc a ) + ^U a ] 



(39) 



where u^j+a (u$ *+ a ) are Hartree-Fock parameters to be 
determined self-consistently. We are interested in real 
parameters which do not depend on the lattice sites, but 
depend on the space directions Ui^+ a = Ui y i+ a = u a . 
Then in momentum space representation, the Hamilto- 
nian (1341) has the form 



H H -F = yE u a^ a -2AMiV + ^£ fc 0+0 fe , (35) 

a fc 

where N is the number of lattice sites and e k is the dis- 
persion of the </>fc-boson (spinon) 

£fc = A — J x u x cos k x — JyUy cos k y — J z u z cos k z (36) 

The free energy of the theory with Hamiltonian (|35[) is 

2T V- 



F =^E u ^ a -2AA/ + ^^ln 



(37) 



where T is the temperature. The equations for the pa- 
rameters w a and A are given by: 



dF 
du a 



= 



dF 



= 



(38) 



To ensure correct definition of the Bose theory (|35|) . 
i.e. to have e k > when the wave vector k runs 
over the first Brillouin zone of a cubic lattice, we have 
to make some assumptions for the parameter A. For 



and the theory is well defined for positive constants u a > 
and > 0. 

We find the parameters u a and /z by solving the sys- 
tem (j3"5)) . For high enough temperatures both fx(T) and 
u a (T) are positive, and the excitation is gapped. De- 
creasing the temperature leads to decrease of /Lt(T). At 
temperature Tq it becomes equal to zero [i(Tc] = 0, and 
long-range excitation emerges in the spectrum. There- 
fore the temperature at which [i reaches zero is the Curie 
temperature. We set — in (f3"5|) and obtain a system 
of equations for the Curie temperature Tq and u a 



ivE ^ 



cos k„i 



k e 



E "Uap a (l — COS fe a ) 



M 



iy 1 



— cos fc a 



(40) 



k e 



- 1 



The results for the Curie temperature Tq as a function 
of doping x are plotted in Figure [5] for different values of 
the Hund's coupling in the absence of Jahn- Teller dis- 
tortions and antifcrromagnctism. One can see that with 
the increase of Ju/t the Curie temperature also increases 
and the maximum of the curves is shifted to lower values 
of x. The ends of the curves correspond to the transition 
from ferromagnetism to antifcrromagnctism and for high 
enough value of J^/t we have ferromagnetism across all 
values of x. This behavior is closely related to the be- 
havior of the spin-stiffness curves (see figure [3]). 
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Figure 5. (Color online) Curie temperature Tc as a function 
of hope doping x for different values of Hund's coupling Jn/t, 
corresponding to the cases depicted in Fig. [3] 



IV. CRITICAL TEMPERATURE IN THE 
A-TYPE ANTIFERROMAGNETIC REGIME 



Let us now turn to calculating the critical tempera- 
ture in the A-type antiferromanetic phase. This phase is 
characterized by two positive and one negative effective 
exchange constants, namely J x > 0, J z > and J y < 0. 
The Hamiltonian then reads: 



i 

-J z Mi-M i+es ] (41) 

We rewrite the vectors in terms of Schwinger 
bosons and use relation (ptU|) for the x and z component, 
while for the y component we use: 

Mi • Mj = -\ (0+0+ - 0+0+) (0.a0 j2 - Mn) 



+ 2 (4*<fe) (0^0^J (42) 

The second term here is a constant and we omit it. 
Introducing a term to ensure the Schwinger bosons con- 
straint, we rewrite the Hamiltonian as 

H = J2 [ Jx {vt+e x v<Vi°') + 

i 

\ J y\ [ ( Ptl l Pt+e y 2 - fiifi+e^l) (<Pil<Pi+e y 2 ~ Vi2fi+e y l) + 

aE(^^-- 2s ) ( 43 ) 

i 

As with the ferromagnetic phase, we treat the four- 
boson interaction in Hartree-Fock approximation, with 
the effective Hamiltonian reading: 



H 



HF 



-\Jy\((vti 



l Pi+e y 2 



2 Yl Jx {VivVi+e*") ( < Pi+e x a' ( P^') + I J v\ (<Pil<Pi+e v 2 ~ V^fi+^l) (<Pil<Pi+e y 2 - <Pi3<fi+e y l) 
i L 

^t+eyl) (<Pim+e y 2 ~ Vi2<Pi+e y \) + {Vil<Pi+e y 2 ~ Vi2Pi+e y l) (pflVi+e^ ~ fa'Pi+eyl)) ~ 

J z ( (Via'Vi+eza') <Pi+e z<T <Pi<T + Wt+e^'Via') vtVi+e^) + A^ (vt^icr - 2s) (44) 



The Hartree-Fock parameters are given by 
{ViaVi+w) = (vt+e^vJ) = u n M = x,z (45) 



■ {}PiWi+e y 2 ~ (Pi2<Pi+e y l) = U y (46) 



I 

the lattice direction a. We then split the Hamiltonian 
into classical and quantum parts 



i?HF — h c i + h q , 
where the classical part is given by 



(47) 



h cl = ±v»N + l -^u 2 y N + ^-u 2 z N - 2XsN, (48) 



and we have again chosen them to be real parameters 

which do not depend on the lattice site i, but depend on while for h q we have 
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i 

J z 



J 



2 x 2 y 2 



(55) 



We can now construct a system of four equations for 



y u z {vt +eza Vio + ^i-^+e^) + XY^PtaPiv ( 49 ) the parameters u s , u y , u z and A: 



This Hamiltonian can be rewritten as 

K = [ £ fe^fea'/'fe CT + iwlk^tk + lkflk<P2k] (50) 



where we have introduced 



Sk = A — J^Ua; cos fc x — J z u z COS k z 



and 



jk = -iJyUy sinky. 
We can easily diagonalize h q to 



with 



(51) 



(53) 



dF dF 
7^ = ° 7n= ' 

OU a OA 

In explicit form the system reads: 

1 ^ £kCOsk x 
mI^ ET (l + 2n fc ) 



Lg^sin^ 



1 ^ £fcCOsfc z 

2s + l = iyfi(l + 2n fc ) 
(52) where is the Bose occupation number 

1 



3 E k /T _ i 



and Ek is the dispersion 



(56) 



(57) 



(58) 



E k = V (A - J x u x cosk x - J z u z cosk z ) 2 - \~f k \ 2 (59) 

The critical temperature is obtained for A = J x u x + 
| Jy\u y + J 2 m z . Substituting into (foTj) we obtain 



£ fe = ^Je 2 -\"f k \ 2 El = |7fe| 2 -£fc (54) £ fe = -4^(1 -cos fc 2 ) + J z u z (1 - cos fc a ) + |J y |%. (60) 

The free energy of the system is then given by 



The dispersion at the critical temperature T = 
then has the form 



E k = J (J x u x (1 - cosk x ) + J z u z (1 - cosfc z ) + \ Jy\uy) 2 - J 2 u 2 sin 2 k v 



(61) 



and E ki = has two solutions for k\ = (0, 0) and k2 = (0, — 0J respectively. Near the zero points, the dispersion 
adopts the form: 



E k ^k z ~ y J x u x \J y \u y kl + J z u z \Jy\u y k 2 z + J 2 ul(^k y T ^ 



(62) 



Thus, the magnons in the antifcrromagnctic phase 
have dispersion that behaves as ~ |fc| at small impulses. 

Using the systems of equations for the Curie and Neel 



temperatures, (|4T))) and (|57[) respectively, we can build 
the phase portraits for the x < 0.5 region. 
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V. RESULTS 

In this section we will apply all the information pre- 
sented so far to construct a realistic phase portrait of 
Lai_ x Ca x MnC"3. As discussed in the introduction, we 
have to describe four phases, namely G-type antifcrro- 
magnet, C-type antifcrromagnet, ferromagnet, and A- 
type antiferromagnet in order of decreasing x. How- 
ever, since the most important transport effects, metal- 
insulator transition and colossal magnetoresistancc effect 
are observed in the ferromagnetic part of the manganitcs 
phase diagram, we will focus only on the x < 0.5 region. 

Before introducing our results, let us first discuss the 
criteria we have used to construct the curves for the crit- 
ical temperatures. Our goal is to be in both qualita- 
tive and quantitative agreement with the experimental 
results, and so we have aimed to reproduce curves with 
the following characteristics: a) the maximum of the fer- 
romagnetic part is observed around x = 0.38; b) the 
maximal value of Tc is around 265 K; c) left and right of 
the maximum the curves are as close to the experimental 
ones as possible; d) transition to A-type antiferromag- 
netic phase occurs at x = 0.08. 

To achieve this, our first step is to determine appropri- 
ate value for Jn/t. As we discussed in section II, values 
lower than Jn/t = 10 cannot reproduce the observed 
physics. Since we also have to consider antiferromag- 
netic term and Jahn- Teller distortions, both of which 
suppress the ferromagnetic phase, we have chosen to work 
with two different values of Jn/i, namely Jn/t = 15 and 
Jn/i = 20. This in turn allows us to establish a lower 
limit upon the value of the hopping parameter t, needed 
to reproduce the observed critical temperatures. For 
J H /t = 15 this value is t = 0.215 cV, and for J H /t = 20 
it is t = 0.20 eV. 

Once we have selected JhA, we turn to the other pa- 
rameters. We have to consider the antiferromagnetic ex- 
change constants J^ F (= Jj(p) ano - ^af ano - the Jahn- 
Tcllcr energy Ejt = g /{2k). For the latter, it is gener- 
ally agreed that its effect decreases with increase of the 
doping level x. For this reason we choose to work with 
distortion, which splits the exchange constants J x and J v 
at x = 0.20 and leads to the phase transition at x = 0.08 
(see fig.E]). This is important, since if the distortion is 
not included, we cannot explain the existence of the A- 
type antiferromagnetic phase. In addition, the absence 
of Jahn- Teller distortion above x = 0.20 means that for 
fixed values of J-n/t the slope of Tc[x] curve is controlled 
by Jl F and J| F . 

To better explain the impact of the antiferromagnetic 
constants on the Curie temperature curves, we have de- 
picted in Figure |6] the phase portraits for three different 
sets of parameters and fixed Jn/t = 15. The blue^ 



cV, and J^ F = 0.000316 eV. The black line corresponds 
0.0140, and cV J% F 



(dashed), gree: 



,26 ; 27 



(dash-dotted) and red^li (dotted) 



lines correspond to the experimental results, while the 
solid lines are obtained using the theoretical calculations 
presented here. The magenta line corresponds to the fol- 



0.000312 

= 0.37 cV, J AF = 0.0142 
0.000421 eV. With increase of J^ F , the 



to t = 0.365 eV, J Z AF = 
cV, and the cyan one to t 
eV, and J£ F 

right end of the curve is shifted to lower values of Tc 
and the maximum of the curve approaches the experi- 
mental value. Decreasing of Jj[ F has the opposite effect, 
and if the value drops below 0.006 cV, the maximum of 
the curve is no longer in the x < 0.5 region. Therefore, 
larger values of Jj[ F are in better agreement with the 
experiment. 

After we have chosen a value of Jj[ F , which ensures that 
the right end and the maximum of our curves are close to 
the experimental ones, we examine the effect of J\ v - Its 
value is important, since, together with the Jahn- Teller 
distortion, it controls the point at which p y reaches zero 
and the phase transition to A-type antiferromagnetism 
occurs. In order to ensure this transition happens at 
x = 0.08, we have to either fix the value of J\ Y and de- 
termine the needed Ejt, or vice versa. With fixed value 
of Jj[ F , ^af a ^ so controls the slope of the curve left of the 
maximum, up to the point of phase transition. Increas- 
ing the value of J\ F lowers the effective constant J y and 
in turn lowers the value of the critical temperature, with 
the effect increasing when we approach x = 0.08. 

The right end on the curves in Figure [6] however is 
not very close to the experimental ones. To bring the 
curves closer, one should further increase the value of 
Jj[ F . However, high values of the antiferromagnetic con- 
stants suppress fcrromagnctism, so in order to further 
increase J AF , we should first increase the value of Jn/i 
to make sure fcrromagnctism persists up to x = 0.08. 

Following the procedure described in the beginning of 
this section, we have obtained the curves for the critical 
temperatures for JhA = 20 (Figure [7]). The larger value 



0.8 n = 1-x 0.7 




Schiffer (exp) 
Cheong (exp) 
Fujishiro (exp) 
T[K] (theory 1 ) 
T[K] (theory 2) 
T[K] (theory 3) 



lowing set of parameters: t = 0.37 eV, J 



AF 



0.00632 



Figure 6. (Color online) Critical temperature as a function 
of hole doping for Jn/t = 15 and: t = 0.365 eV, J A f = 
0.0140, eV Jaf = 0.000312 eV (black line); t = 0.37 eV, 
Jap = 0.0142 eV, J1f = 0.000421 eV (cyan line); t = 0.37 
eV, Jlp = 0.00632 eV, Jaf = 0.000316 eV (magenta line). 
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for the Hund's constant allows us to work with larger val- 
ues for the antiferromagnetic ones, namely J' 



AF 



0.0206 



eV and J^ F 



0.00176 eV, which in turn results in curves 
that are in better agreement with the experiment. To 
better illustrate the effect of the last remaining parame- 
ter, the Jahn- Teller distortion, we have used the same set 
of parameters for all the curves. Thus they are identical 
up to the point where distortion effects set off (around 
x = 0.2). As in figure [6l blue (dashed), green (dash- 
dotted) and red (dotted) curves correspond to the exper- 
imental results. 

We have worked with three different types of distor- 
tion, which are shown in the inset of Figure [7] Their 
crossing point represents the value of Q{n), which we 
need in order to have transition to A-type antifcrromag- 
nctism at x = 0.08 (i.e. to have J y — 0). This value 
is determined by the choice of Jn/t and J\p, so m this 
case it is equal for all three curves. In Figure [7J the black 
curve is a reference curve that corresponds to distortion, 
which increases linearly from x = 0.2 to x = 0. The ma- 
genta curve corresponds to distortion Q(n), obtained by 
minimizing the fermion part of the free energy^. Both 
the black and the magenta curves however, while giving 
values for Tpf very close to the experimental ones, de- 
crease when x approaches zero. To have the same behav- 
ior as the experimentally observed curve, the distortion 
has to grow rapidly for small values of x. Such behav- 
ior might be explained if we consider anharmonic terms 
in the phonon Hamiltonian. The cyan curve in Figure [7J 
corresponds to phcnomenologically fitted distortion, such 
that the resulting curve for Tjv has the same behavior as 
the experimentally observed one. It only deviates from 
the magenta curve in the A-type antiferromagnetic phase. 



0.8 n = 1-X 0.7 




Schiffer (exp) 
Cheong (exp) 
Fujishiro (exp) 
T[K] (theory 1) 
T[K] (theory 2) 
T[K] (theory 3) 



Figure 7. (Color online) Critical temperature as a function 
of hole doping for Jn/t = 20, t = 0.462 eV, J AF = 0.0206 
eV, Jaf = 0.00176 eV, and different types of Jahn- Teller dis- 
tortion. Inset: Jahn- Teller distortion as a function of carrier 
density (hoping). 

One can see that near the critical doping value, our 
curves start to rapidly decrease. The reason behind this 
effect is the decreasing value of J y . As its value ap- 



proaches zero, we are in effect describing the 2D case. 
Since our method of calculation is in agreement with 
the Mermin- Wagner theorem^, we correctly obtain zero 
temperature at the critical doping value x = 0.08. How- 
ever this is not the observed experimental behavior and 
is a result of the limitations of our method, which can be 
avoided if we consider next to nearest neighbor correc- 
tions to the effective Hamiltonian (|26[) . For this reason, 
in our phase portraits we have shaded the region around 
x = 0.08, where nearest neighbor contributions result 
in diminishing value of J y . If one accounts for next to 
nearest neighbor corrections, the curves will continue to 
smoothly decrease past the critical point, which is repre- 
sented by the dashed continuation lines near the critical 
point in Figure [7] 



VI. SUMMARY 

Starting from the well known Double Exchange Model 
and supplementing it with antiferromagnetic and Jahn- 
Teller terms, we have derived effective Heisenberg-type 
model for a vector, which describes the local orienta- 
tion of the total magnetization. We have then used 
Schwingcr-bosons mean-field theory in Hartrcc-Fock ap- 
proximation to calculate the critical temperatures in the 
ferromagnetic and antiferromagnetic regimes. This tech- 
nique of calculation is in agreement with the Mermin- 
Wagner theorem. We have then shown that the combi- 
nation of these two ingredients provides results for the 
critical temperatures, which are in very good agreement 
with the experimental results. 

We have argued that, in order to explain all the ob- 
served phases, one has to consider values for the Hund's 
coupling as large as Jji/t = 20. Another key point is 
to use antiferromagnetic constants, which depend on the 
lattice direction. Indeed, best agreement with the exper- 
imental results is observed when the value of </j[ F is much 
larger than J^ F and J^ F - Based on the agreement with 
the experiment, in Table U we provide estimation on the 
values of the model parameters, which best describe the 
observed physics. We have summarized the "best fit" val- 
ues for both cases we have examined, namely Jn/t = 15 
and J s /t = 20. 

While the method we have used is in very good agree- 
ment with the experimental results, it can still be im- 
proved. The inclusion of next to nearest neighbor correc- 
tions is needed to avoid the rapid decrease of the critical 
temperature curves near the point of phase transition. 



Table I. Estimation for the model parameters in eV for 
Jn/t = 15 and Js/t = 20. 







t 


•Jaf 


V 

AF 




Jn/t = 


15 


0.37 


0.000421 


0.000421 


0.0142 


Jn/t = 


20 


0.462 


0.00176 


0.00176 


0.0206 
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Anharmonic terms in the phonon Hamiltonian should portance of the tolerance factor and related structural 
also be considered and will possibly provide better agree- details such as Mn — O — Mn bond is another thing we 
ment with the experiment in the low x limit. The im- have not considered. 
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